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Abstract. We present the first measurement of reconstructed jets in ultra-relativistic 
heavy ion collisions. Utilizing the large coverage of the STAR Time Projection 
Chamber and Electromagnetic Calorimeter, we apply several modern jet recon- 
struction algorithms and background subtraction techniques and explore their 
systematic uncertainties in heavy ion events. The differential spectrum for inclu- 
sive jet production in central Au-|-Au collisions at ^/sNN ~ 200 GeV is presented. 
In order to assess the jet reconstruction biases, this spectrum is compared with 
the jet cross section measured in ^/s — 200 GeV p-f p collisions scaled by the 
number of binary N-N collisions to account for nuclear geometric effects. 

1 Introduction 

Strong suppression of inclusive hadron distributions and di-hadron correlations at high pT is 
observed at the Relativistic Heavy Ion Collider |l|2j . This provides direct evidence for partonic 
energy loss in dense matter generated in ultra-relativistic heavy ion collisions. However, such 
measurements suffer well-known geometric biases due to the competition of energy loss and 
fragmentation - the leading particle spectrum is dominated by relatively low energy jets that 
happen to lose little energy in the medium and fragment into higher px particles [3]. These 
biases can be removed if the underlying partonic kinematics are reconstructed in an unbiased 
way, independent of the fragmentation details - quenched or unquenched. An unbiased jet 
reconstruction measurement in heavy ion collisions would give access to the full spectrum 
of fragmentation topologies without geometric biases, enabling full exploration of quenching 
dynamics. In addition, fully reconstructed jets allow the measurement of qualitatively new 
observables such as jet shapes, fragmentation functions, and energy flow. 

In this article, we present the first systematic measurement of direct jet reconstruction in 
ultra relativistic heavy ion collisions. We study the effects of high multiplicity underlying events 
of reconstructed jets by applying several modern jet reconstruction algorithms and background 
subtraction techniques that are appropriate to heavy ion collisions. See [4| for the accompanying 
discussion about jet fragmentation studies in heavy ion collisions also presented for the first 
time during the Hard Probes 2008 meeting. 



2 Analysis and Techniques 

2.1 Event Selection and Terminology 

This analysis utilizes Au-|-Au data at ^/snn — 200 GeV recorded by STAR for the minimum 
bias and enhnanced high pT collisions during the year 7 RHIC run. The most central (0-10%) 
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Au+Au collisions are selected, using the multiplicity measurement performed with the STAR 
Time Projection Chamber (TPC). Two event sets were analyzed, based on fast on-line trigger 
configurations: 

(i) MB-Trig, (minimum-bias trigger) utilizing the coincidence between the two calorimeters at 
beam rapidity (Zero Degree Calorimeters) with signals in each greater than ^ 40% of the most 
probable amplitude for a single neutron, and 

(ii) HT-Trig (high-tower trigger) that satisfies the MB-Trig conditions and the additional re- 
quirement of a 2-tower EMC clusters having at least a 7.5 GeV energy deposition. 

For this study we use 3 million MB-Trig events, corresponding to 300 thousand 0-10% most 
central events. A total of 80 million MB-Trig events were recorded by STAR during the 2007 
Au + Au run, but at present only the event set we use have been fully reconstructed off-line. 
The HT-Trig is designed to enhance the recorded rate of high px photons and electrons. It 
may also serve to enhance the recorded rate of jets, though potentially with a bias that we 
investigate in this study. The HT-Trig dataset, which samples the full integrated luminosity of 
Run 7, has been reconstructed in its entirety, corresponding to 500 /ib^^. 

In order to assess jet reconstruction energy resolution, background subtraction, efficiency 
and acceptance, we utilize Monte-Carlo model studies based on Pythia 8.107 [5]. Pythia events 
with hight Et jets are generated in three different configurations: 

(i) PyTrue: Pythia isolated jets including all particles except neutrinos. Jets are reconstructed 
using the PYTHIA internal jet finder, PyCell, for the cone algorithm, and FastJet for the 
sequential recombination algorithm. 

(ii) PyDet: Pythia isolated jets (parameterized detector response level) reconstructed using the 
algorithms we apply to real data. 

(iii) PyEmbed: PyDet that are embedded in a background of real Au-|-Au 0-10% central events, 
with jets reconstructed using the algorithms we apply to real data. 

We use several jet reconstruction algorithms to study the systematics. See next section for 
the detailed discussion of the algorithms. In all of these algorithms, only the highest energy jet 
per event is selected as the reconstructed jet. 



2.2 Jet Reconstruction Algorithms and Background Corrections 

Jets are reconstructed by combining charged particles from the STAR TPC and the neutral 
energy from the Electromagnetic Calorimeter (EMC). TPC and EMC detectors cover full az- 
imuth {0 < (j) < 27r) and mid rapidity (—1 < 77 < 1) of the events. Corrections are applied for 
double-counting of energy due to hadronic energy deposition in the EMC and to electrons. Two 
kinds of jet reconstruction algorithms are utillized; seeded cone (leading order high seed cone 
(LOHSC)) and sequential recombination {Kt and Cambridge/ Aachen). For an overview of jet 
algorithms in high energy collisions, see [6|7|8|9|10] and references therein. In the following, we 
briefly discuss the two kinds of algorithms and the corresponding underlying event subtraction 
that we use. 

The main method of the cone algorithm is to combine particles in ?/ — </) space with their 
neighbors within a cone of radius R (i? = y^A(j)^ + Arf) . To optimize the search and effective- 
ness of jet finding in p + p and e+ + collisions, splitting, merging, and iteration steps are 
used. In order to avoid instabilities in cone-finding due to large heavy ion background, we use 
a simple seeded cone without iteration or split-merging steps, with cone radius R = 0.4 and 
minimum seed of 4.6 GeV. A relatively small cone size is chosen to suppress the underlying 
heavy ion background |11I12| . In p + p collisions ~ 80% of the jet energy is observed to be 
within R~ 0.3 for 50 GeV jets in the Tevatron [13]. However, broadening of the jet fragmen- 
tation due to quenching in the medium formed in heavy ion collisions may reduce the fraction 
of the measured energy in a given cone size and needs further exploration. To further reduce 
the heavy ion background, the minimum accepted transverse momentum of charged particles 
and the transverse energy of the calorimeter cells (p™*) is varied between 0.1 to 2 GeV. This 
threshold cut does not remove all the background contamination on the jet energy and addi- 
tional subtraction is needed. The residual background is corrected based on the out-of-cone 
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energy for the same averaged over the STAR acceptance but measured on an event-wise 
basis, and scaled to the cone area. 

The sequential recombination algorithms have been used extensively in the Tevatron as they 
are coUinear and infrared safe [13 14 . In these types of algorithms, arbitrary shaped jets are 
allowed to follow the energy flow resulting in less bias on the reconstructed jet shape than with 
cone algorithms [8]. The Fast Jet package for sequential recombination algorithms was used 
for this analysis [15116] . The Fast Jet package includes kr, Cambridge/ Aachen, anti-fcr, and 
an interface to external jet finders such as SisCone via a plugin mechanism. For infrared and 
collinear safe algorithms an active area (Aj) of each jet is estimated by filling an event with 
many very soft particles and then counting how many are clustered into a given jet. If the 
underlying event is distributed uniformly in rj and (f> then this noise density can be subtracted 
from the measured jet energy on an event-by-event basis to correct for the background energy 
underlying the jet. In simulations, this correction is observed to recover most of the jet's energy 
when they are reconstructed in pile up and heavy ion backgrounds [8]. 
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Fig. 1. Jet area and fluctuations are reconstructed utilizing the Fast Jet algorithm [8] for real jets 
in minimum bias triggered 0-10% central Au + Au collision (MB-trig), in Pythia isolated jet events 
embedded in real central Au + Au events (PyEmbed) and in Pythia isolated jet events (PyTrue). 



Figure [T] shows the jet area and the fluctuations obtained directly from Fast Jet framework 
with the algorithm for the PyTrue, PyEmbed and MB- Trig jets. See the insets of the 
histograms for the applied parameters. Pythia jets embedded in real Au-|-Au background events 
are observed to have the same area and fluctuations as jets from real Au-|-Au events. This can 
be seen from the comparison of solid blue and dashed red histograms in Figured] A reduction 
in the jet area and an increase in fluctuations are observed for the MB-Trig jets when compared 
with the PyTrue jets. This reduction of the jet area relative to PyTrue is well understood for 
sequential recombination algorithms on theoretical grounds 8J. 

The inclusive jet spectrum from p + p collisions at y/s — 200 GeV {p + p STAR) is used 
to make comparisons with the jet results in heavy ion collisions. Jets in p + p collisions are 
measured the same way as in Au + Au collisions, utilizing the STAR TPC and EMC and 
correcting for missing and double counted energy [17j . However for the p + p case, a mid-point 
jet cone algorithm with splitting and merging steps is used. The inclusive jet spectrum agrees 
well with the Next-to-leading order perturbative QCD calculation [20] . 
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2.3 Energy Resolution 



The energy resolution for jet reconstruction with various algorithms has been studied with 
isolated jets simulated with Pythia [5 18 19J. Figure [5] shows the event by event comparison of 
PyTrue, PyDet and PyEmbed from LOHSC, kr and Cambridge/ Aachen algorithms. Applied 
cuts, jet energy and the algorithm used, are specified in the figures. For all three algorithms, 
a shift of median due to un-measured particles, primarily neutrons and K^, and the applied 
Pt cut (hence loss of jet energy), is observed for the AE = EpyDct — EpyXruo histograms. 
Background effects are simulated using Pythia jets that are embedded in real Au + Au events. 
The distributions in Figure [5] is convoluted with the true jet spectrum to produce the observed 
jet spectrum. The effect of the heavy ion background on the jet energy can be seen in the 
AE = EpyEmbod — EpyTruo distributions. The algorithms show slightly different behavior: the 
LOHCS algorithm has more entries with AE > while the kx algorithm shows fewer entries at 
AE > 0, but has average AE < , i.e. a overall oversubtraction of the background due to the 
back-reaction of the background (smaller average area, see Figure [1] and [S]). A positive AE in 
this distribution can distort the measured inclusive jet spectrum substantially, increasing the 
apparent yield at high Et and resulting into a harder spectrum. A correction to the spectrum 
must be applied to account for this effect. 




LOHSC 


1600 
1400 


; R=c 


Et=35+5 GeV 


1200 




larv 


1000 
800 




STAR 
Prelimir 


eoo 

400 
200 




D ID 20 3 

AE 




^0 



=1 GeV 



CAMB 

Et=35±5 GeV 



Fig. 2. Distributions showing energy resolution; black is AE 

EpyEmbcd — EpyTruc and green is AE — EpyEmbcd ^ EpyOet. 
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The influence of energy resolution on jet spectrum can be observed with PyDet, PyEmbed 
and PyTrue distributions shown in Figure [3] for the LOHSC algorithm for p™* of 0.1, 1 and 
2 GeV. In the smallest p™* case, a large difference between PyEmbed and PyDet is observed, 
while in the p™* = 2 GeV case when the background contributions is the smallest, PyEmbed 
is similar to PyDet. The kick in the jet spectrum for PyEmbed is observed in the p™* = 0.1 
GeV, as expected from the tail at positive AE due to large background in red distributions in 
Figure [21 When the px threshold is increased, the background fluctuations are reduced and the 
enhancement in the spectrum relative to the case without background is reduced to a negligible 
level. With this larger p^?** = 2 GeV threshold requirement, jet reconstruction in 0-10% most 
central Au-f-Au collisions is similar to that of p+p collisions. However, while a reduction in the 
background fluctuations is achieved, a reduction in the measured jet energy (p™* dependent 
bias) is introduced. Similar effects, though smaller in magnitude, are also observed for the Kx 
and Cambridge/ Aachen algorithms. 



2.4 Jet Spectrum and Corrections 



The correction factors of the jet spectrum are estimated using PyEmbed. The comparison of the 
spectrum of the jets of PyTrue and PyEmbed are shown in the left panel of Figured The Ex 
dependent ratio of PyEmbed to PyTrue is shown in the right panel of Figure ID A polynomial 
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Fig. 3. Inclusive jet spectrum for PyDet, PyEmbed and PyTrue using the LOHSC algorithm with 
three different pr cuts on track momentum and calorimeter cell energy. Note the lower threshold on 
generated jet energy e!^^'^^^" > 5 GeV, which affects the reconstructed spectrum up to Et = 20 GeV. 



function fit to the distribution, shown as the dashed histogram, is used as a muhiphcative 
correction to the inclusive spectrum. Table[T]shows the inclusive jet spectrum correction factors 
for various p™* values. For the sequential clustering algorithms, the correction factors are closer 
to unity. 
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Fig. 4. The inclusive jet spectra from PyTrue and PyEmbed constructed with the LOHSC are shown 
on the left and their ratio is shown on the right. 



The corrected inclusive jet spectrum for the LOHSC algorithm for the p™* = 1 GeV is 
presented in Figure O The p™* = 1 GeV is selected as it corresponds to correction factors close 
to unity. At this p™* the competing effects of energy loss due to momentum threshold cut and 
the kick up in the jet spectrum due to the positive tail of energy resolution cancel each other. 
The filled triangles are for the MB-Trig data set and are corrected for resolution, acceptance 
and efficiency. 
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Table 1. Correction factors for the inclusive jet spectrum, for different reconstruction algorithms and 
values of p™*. The range of values given indicates the correction factor variation from the lowest to 
highest jet Et shown in the figures. 
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2.5 Comparison to p + p Spectrum 

In order to assess the biases in the jet spectrum reconstructed in central Au-|-Au collisions, 
we compare to the spectrmxi measured in p-l-p collisions. To account for nuclear geometric 
effects we scale the p+p spectrum by Nsinary, the number of binary nucleon+nucleon collisions 
equivalent to a central Au-|-Au collisions, as calculated by a Glauber model ^J. The cross 
section for hard processes is expected to scale with Neinary if no nuclear effects are present. 
In the case of jet reconstruction, Neinary scaling is expected if the reconstruction is unbiased, 
i.e. the jet energy is recovered independent of the particular mode of fragmentation, even in 
the presence of strong jet quenching. The NBinary scaled jet spectrum from p + p collisions 
is shown in filled squares [17]. The yellow band represents the systematic uncertainity of the 
p + p jet measurement. Heavy ion jet spectrum is observed to agree with NBinary scaled p+p 
measurement within the ~ 50% systematic uncertainty of the normalization. 

Open circles in Figure[5]show uncorrected jet spectrum from HT-Trig data which is substan- 
tially lower then the corrected MB-Trig spectrum. The correction factors of energy resolution, 
efficiency and acceptance for the HT-Trig are expected to be small. A large trigger bias due 
to the additional 7.5 GeV energy deposition in EMCAL requirement fot HT-Trig relative to 
MB-Trig is seen to persist to at least to 30 GeV. Further statistics of MB-Trig data is needed 
to assess the bias at high px- 

Figure [6] shows the comparison of inclusive jet spectra from the MB-Trig Au + Au data and 
the Nginary scaled p+p for the variations of px threshold cut for LOHSC and kx algorithms. 
As can be seen from this Figure, the agreement between NBinary scaled p+p and MB-Trig 
measurement is good for the lowest value of pt cut. It is also seen to be poorer with the larger 
Pt threshold cut. This suggests that p™* introduces biases which are not fully corrected using 
Pythia as the fragmentation model, and may be an indirect indication of modified fragmentation 
due to quenching. 



3 Summary and Outlook 

We have presented the first full reconstruction of jets in 0-10% most central heavy ion collisions 
at RHIC energies. Systematics of the underlying heavy ion background subtraction is studied 
via utilizing various algorithms and consideration of the jet area. The Neinary scaling is observed 
for the least-biased cuts with the given ~ 50% systematic uncertainty of the p + p jet spectrum 
measurement. An unbiased jet reconstruction in heavy ion collisions with MB-Trig data appears 
to be feasible. However, spectrum corrections are currently based on model calculations using 
Pythia fragmentation. This aspect, together with background subtraction techniques, spectrum 
variations due to cuts and reconstruction algorithms, must be investigated further in order to 
assess the systematic uncertainties of this measurement. 

The LHC will begin operations in 2008, with the first heavy ion run expected in late 2009. 
At the LHC the heavy ion background is expected to be larger than at RHIC, but there will be 
copious production of very energetic jets, well above background |22j . The large kinematic reach 
at the LHC may provide sufficient lever-arm to map out the QCD evolution of jet quenching 
[23] . The comparison of full jet measurements in the different physical systems generated at 
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Fig. 5. Jet yield per event vs Et for 0-10% central Au+Au collisions, compared to the distribution from 
p + p collisions scaled by Nsinary [IZ]. Triangle symbols are from MB-Trig and corrected for efficiency, 
acceptance and energy resolution. Open circles are from HT-Trig and not corrected for trigger bias. 
Only statistical error bars are shown for the Au + Au data. Filled black squares are the distribution 
from p + p collisions, scaled by N Binary Band represents the systematic uncertainty of the p + p 
measurement. 



IIJ- 1 


r Au+Au 0-10% 




MB-Trig 




10 ' 

10* 


r V^: ■ Km Scaled p+p 


10^ 
10' 


_ 0-10 % 

^ LOHSC 




10* 


, pTcLt>0.1GeV ^ 




10' 

10* 


r STAR Preliminary 

1 

10 20 30 40 





Au+Au 0-10% 



2*10' 
10" 



10' 



Au+Au 0-10% 

^ MB-Trig 

■ Nbin Scaled p+p 



\\\ 



KT 

pTcut> 0.1 GeV 



10^ 



0-10% ' 

LOHSC 



Au+Au 0-10% 

A MB-TrIg 
■ Nan Scaled p+p 



pT cut > 1 GeV 

STAR Preliminary 



E, [(SeVl 



z'io' 

10' 



r Au+Au 0-10% 


1 V ■ 


MB-Trig 
Nh„ Scaled p+p 








r KT 






r pTcLt> 1 GeV 

■ 








10^ r 



10* 



Au+Au 0-10% 

A MB-Trig 
■ N^in Scaied p+p 




LOHSC 

pT cut > 2 GeV 



STAR Prelim 



40 SO 
E, [GeV) 



Au+Au 0-10% 

A MB-Trig 
'W ■ Nbin Scaled p+p 



KT 

pT cut > 2 GeV 



40 50 

■ - c, Ifevi E, |SeV| 

Fig. 6. The pr threshold dependent comparison of inclusive jet spectra in heavy ion and p-|-p collisions. 
The filled triangle symbols are from MB-Trig and filled squares are from Nsinary scaled p + p collisions. 
Upper panel is for the jet yield per event vs Et reconstructed with LOHSC algorithm and the lower 
panel is for the same distributions with the kr algorithm. In addition to pr cuts shown in the figures, 
seed is 4.6 GeV for the LOHSC and R is set to 0.4 for both algorithms. Only statistical error bars are 
shown. 



RHIC and the LHC will provide unique and crucial insights into our understanding of jet 
quenching and the nature of hot QCD matter. 
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